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Pairing gaps for fermionic atoms in harmonic oscillator traps are calculated for a wide range of 
interaction strengths and particle numbers, and compared to pairing in nuclei. Especially systems 
where the pairing gap exceeds the single level spacing but is smaller than the shell splitting fiu), 
are studied which applies to most trapped Fermi atomic systems as well as to finite nuclei. When 
solving the gap equation for a large trap with such multi-level pairing, one finds that the matrix 
elements between nearby harmonic oscillator levels and the quasi-particle energies lead to a double 
logarithm of the gap, and a pronounced shell structure at magic numbers. It is argued that neutron 
and proton pairing in nuclei belongs to the class of multi-level pairing, that their shell structure 
follows naturally and that the gaps scale as ~ A -1 / 3 - all in qualitative agreement with odd-even 
staggering of nuclear binding energies. Pairing in large systems are related to that in the bulk limit. 
For large nuclei the neutron and proton superfluid gaps approach the asymptotic value in infinite 
nuclear matter: A ~ 1.1 MeV. 



PACS numbers: 03.75.Fi, 21.65.+f, 74.20.Fg, 67.60.-g, 05.30.Fk 
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I. INTRODUCTION 

Pairing in Fermi systems is central for understanding 
superconductivity in condensed matter physics [1], su- 
perfluidity and glitches in neutron stars, excitations spec- 
tra and odd-even staggering of binding energies in nuclei 
[2], metallic clusters [3] and superconducting grains [4]. 
New insight into the general properties of such Fermi sys- 
tems can now be obtained from experiments with atomic 
gases [5] which recently have been cooled down to de- 
generate temperatures around 10% of the Fermi temper- 
ature [6-9]. The interactions between trapped 6 Li and 
40 K atoms have simultaneously been tuned above a Fes- 
hbach resonance, where they become strongly attractive, 
in order to produce optimal conditions for superfluidity 
[10]. Atomic gases have widely tunable number of parti- 
cles, densities, interaction strengths, temperatures, spin 
states, and other parameters which holds great promise 
for a more general understanding of pairing phenomena 
in solids, metallic clusters, grains, nuclei and neutron 
stars. 

In the experiments with 6 Li atoms the anisotropic ex- 
pansion after sudden release from the trap is as predicted 
from hydrodynamics [11]. This is compatible both with 
a strongly interacting superfluid and collisional hydrody- 
namics. The scattering length is attractive a < and 
large such that densities above p ^> l/|a| 3 are achieved 
and found mechanically stable against collapse. The in- 
teraction energy per particle, like Fermi energy, scales as 
p2/3 [6-8] a s predicted theoretically [12,13]. The super- 
fluid gaps in bulk are also expected to be of order the 
Fermi energy. 

In nuclear physics pairing is observed directly in the 
odd-even staggering of binding energies, i.e., that nuclei 
with an even number of protons or neutrons are more 
strongly bound than nuclei with corresponding odd num- 
bers by a pairing gap of order ~ 1 MeV. The same pair- 



ing gap also determines the excitation energies of ground 
state nuclei. On average the pairing energies in nuclei 
decrease with the number of protons and neutrons modu- 
lated with a distinct shell structure, i.e., they are smaller 
near closed shells. It would be most interesting if the 
pairing gaps of fermionic atoms in harmonic oscillator 
traps would display similar shell structure and scaling 
with the number of particles as nuclei. 

The aim of this paper is to calculate pairing gaps 
in ultracold atomic Fermi gases in harmonic oscillator 
traps and in nuclei, which are then compared to data on 
odd-even staggering of nuclear binding energies. For the 
atomic traps existing pairing gap calculations [15,16] are 
extended to systems, where level spectra and shell effects 
play an important role. It is shown that such pairing 
mechanisms are similar in nuclei, which to first approxi- 
mation can be considered as a finite system of fermions 
in a harmonic oscillator field with attractive delta func- 
tion interactions. Another important common feature is 
the anharmonic fields that lead to splitting of the sin- 
gle particle states which reduce the pairing to the levels 
closest to the Fermi surface and leads to a distinct shell 
structure. For a wide range of interaction strengths and 
number of particles the trapped atomic clouds are pre- 
dicted to display similar scaling and shell structure as 
nuclei as also seen in the experimental data on neutron 
and proton pairing. The scaling with particle number 
and the continuum or BCS limits are calculated for both 
large traps and nuclei and thereby the pairing gaps are 
estimated for nuclear matter, which also gives an idea 
about the superfluid gaps in neutron star matter. 

The paper is organized as follows. In section II the 
basic properties of interacting fermions in harmonic os- 
cillator traps are given, in particular the single particle 
level spectra. Pairing is treated in section III with indi- 
vidual subsections devoted to each pairing regime where 
the last three subsection outline the novel bulk, strongly 
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interacting and multi-level pairing regime. In section IV 
we turn to nuclei and show that they belong to the class 
of multi-level pairing. 



II. DILUTE FERMI GASES IN HARMONIC 
TRAPS 

We treat a gas of N fcrmionic atoms of mass to in a 
harmonic oscillator (h.o.) potential at zero temperature 
interacting via a two-body interaction with attractive s- 
wave scattering length a < 0. We shall mainly discuss a 
spherically symmetric trap and a dilute gas (i.e. where 
the density p obeys the condition p|a| 3 <C 1) of particles 
with two spin states with equal population. The Hamil- 
tonian is then given by 



H = 



N 

E 

i=l 
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(1) 



with the effective coupling g = Airh 2 a/m. For a large 
number of particles N at zero temperature the Fermi 
energy is for a non-interacting system 



E F = (n F + 3/2)Jku ~ (37V) 1/3 Hw 



(2) 



where n F ~ (3-/V) 1 / 3 is the h.o. quantum number at the 
Fermi surface. The h.o. shells are highly degenerate with 
states having angular momenta I = rip, rip — 2, ...,1 or 
due to the 17(3) symmetry of the 3D spherically sym- 
metric h.o. potential. However, interactions split this 
degeneracy. In the Thomas-Fermi (TF) approximation 
(see, e.g., [17]) the mean-field potential is 



U(r) = 
the Fermi energy 
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The density can be determined from the Kohn-Sham en- 
ergy density functionals. In a dilute gas the mean field is 
small as compared to the Fermi energy and its effect on 
the density can ignored in the following. Thus 



p(r) = k F (r)/3ir 2 ~ p Q (l - r 2 /R 2 TF ) 



3/2 



(5) 



inside the cloud 



< R 



TF 



a osc ^/2np~+~5, where 

a osc = ^/h/moj is the oscillator length, and po = (2n F + 
3) 3 / 2 /3ir 2 a 3 osc the central density. [18] 

The splitting of the h.o. shells degenerate levels / = 
n Fl n F — 2, 1 or in the shell n F by the mean-field 
potential can be calculated perturbatively in the dilute 
limit. An excellent approximation for the radial h.o. 
wave function with angular momentum I and (n F — I)/ 2 



radial nodes in the h.o. shell when n F ^> 1 is the WKB 
one [14] 



2 sin(kir + 6) 
Jtt k t r 



(6) 



between turning points [14]. The phase 9 will not be 
important in the following. The WKB wave number fc; (r) 
is given by 



kf(r) 



2n F 
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The single particle energies are 
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3ir a 



— (2n F + 3f/ 2 hLo 



U{r)\K nFl {r)\ 2 r 2 dr 
4 1 1(1 + 1) 



37r 47r n F 



(8) 



The latter result is calculated from the overlap between 
the mean field as given by Eqs. (3) and (5) and the 
WKB wave-functions as given by Eqs. (6) and (7) [15]. 
It is exact within WKB for 1 -c n F and it compares well 
to numerical results for all 1. For attractive interactions 
(a < 0) the lowest-lying states have small angular mo- 
mentum, which is opposite to nuclei. 

An important quantity for pairing is the supergap 
which was introduced in Ref. [15] and will be discussed 
below 



_ 32V2njr + 3 \a\ 
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In comparison the total shell splitting in a shell n F is 
fromEq. (8) 



D 



enj^=0| 



5(n F + 3/2) 
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(10) 



The pairing depends crucially on the shell splitting D 
and the supergap G as will be shown in the following 
section. 
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FIG. 1. Diagram displaying the regimes for the various 
pairing mechanisms (see text) at zero temperature in h.o. 
traps vs. the number of particles N = n%/3 and the inter- 
action strength a. The dotted lines indicate the transitions 
between single-shell pairing A = G, multi-level Eq. (30), sin- 
gle-level Eqs. (15-17), and multi-shell pairing Eq. (24). At 
the dashed line determined by 2Gln(7riF) = hui the pairing 
gap is A ~ Tiu), and it marks the transition from multi-shell 
pairing to bulk superfluidity Eq. (25). The pairing gap is 
A = 0.54Z5f above the full line p|a| 3 > 1, which separates the 
dilute from the dense gas. 



III. PAIRING 

Pairing in small systems as nuclei [16], metal clusters 
[3] and superconducting grains [4] appears in odd-even 
staggering of binding energies, i.e. that systems with 
even number of particles are more strongly bound than 
odd ones. In nuclei the pairing is also responsible for 
superfluid effects in collective motion [19]. 

Pairing in finite systems is described by the 
Bogoliubov-de Gennes equations [21] and take place be- 
tween time-reversed states. As shown in [16] these states 
can be approximated by h.o. wave functions in dilute h.o. 
traps as long as the gap does not exceed the oscillator en- 
ergy, A^,hu>. The gap equation at zero temperature then 
reduces to 



21' + 1 A n /j, : . 2 , 



Tj, 4vr E n , v J 



with quasi-particle energies 

E n n> = \J {e n ' i' - M) 2 + A 2 ,,, ■ 



(11) 



(12) 



The cut-off n^flnp in the sum of the gap equation models 
as a first approximation the more rigorous regularization 
procedure described in Ref. [22] that is required for a 
delta-function pseudo-potential. 

The chemical potential [i is related to the number of 
particles by 



N = 2^2(21' + !) 



e n 'l' 



E„ 



(13) 



For weak interactions the chemical potential lies within 
A n i from the level energy at the Fermi surface e n i , except 
for closed levels where the chemical potential lies between 
the closed and next open level. For closed shells N = 
(«F + l)(np + 2)(np + 3)/3, the gap vanishes for inter- 
actions below a critical value, G < G c — Tiu /2\n{^np). 
[16] 

The gap at the Fermi surface is most important and 
it is generally proportional to the critical temperature as 



discussed in the appendix. For half filled level (or shell), 
where /i = e nF i, the gap is equal to the quasi-particle 
energy E np i — A np j. We will in the following refer to 
this gap as A = E npl = A„ pi . 

The h.o. traps provide a physical system in which 
the connection between these methods of pairing calcula- 
tions can be demonstrated and analytical results can be 
given in various limits of interaction strength and particle 
number. The associated pairing regimes are displayed in 
Fig. 1 and will be discussed in the following subsections. 
The first three have been discussed in earlier publications 
[15,16] but contains additional details of derivations and 
sets the notation. In the remaining three subsections, 
which describe new pairing regimes, pairing gaps are cal- 
culated and the transitions to the other regimes are dis- 
cussed. The last subsection treat the important regime 
of multi-level pairing which will be shown to be especially 
relevant for nuclei in the following section. 



A. Single-level pairing 

For very weak interactions pairing take place only be- 
tween time-reversed states V , /m( r i)^;-m( r 2) within the 
Mevel at the Fermi surface, i.e. I' = I. Maximal pairing 
between two particles in the level is achieved with the 
Cooper pair wave function 

0o(n F , = "^(Iml - m|00)V* m (ri)V*-m(r 2 ) • (14) 

The pairing energy between only two particles in the level 
is twice the quasi-particle energy, and so 



E, 
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drr 2 nt Fl {r) 



(15) 



It follows from the seniority scheme (see, e.g., [23]) that 
this is generally the quasi-particle energy and excitation 
energy for any even number of particles in the level. Solv- 
ing the gap Eq. (11) also gives E np i = A for any number 
of particles in the level. 

The pairing energies were given in Ref. [15]. For the 
top levels I ~ np the h.o. wave functions are sufficiently 
simple that the pairing energy can be calculated exactly. 
For example, n 2 npl=nF = (2/T(l + 3/2))r 2 ' exp(-r 2 ), (in 
units where a osc = 1), which inserted in Eq. (15) gives 
the gap 



E nf 



-huj . 



(16) 



The h.o. wave function for I = can be approximated 
by the WKB one of Eq. (6) when np ^> 1. Inserting in 
Eq.(15) gives the pair gap [14] 



E nF l=Q — 



V2 



a 



irwnp a osc 



(17) 
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Around midshell, i.e. away from closed (1 <C I) and 
open (1 -C uf — I) shells, the WKB wave functions must 
be replaced near their turning points by Airy functions 
which removes a logarithmic singularity and in turn leads 
to a logarithm ln(Z). We find [24] 



E nF i = ^<<&o(nF)|<5 3 (ri - r 2 )\%(n F )} 



■n A Jn^ a osc 



l«:Z<n F , (18) 



to leading logarithmic accuracy. 

As shown in Fig. 1 the single-level pairing regime ap- 
pears for sufficiently weak interactions and large number 
of particles (rip £50) such that D 3> G. The correspond- 
ing pairing gaps can be seen in Fig. 2 at large np. 



B. Single-shell pairing 

When rip is sufficiently small the mean-field splitting 
D is smaller than the pairing gaps and the pairing acts 
between all states in an oscillator shell and not just in a 
single I multiplet, i.e. the full SU(3) symmetry is effec- 
tively restored as compared to the SU (2) symmetry of 
the single I multiplet. This enhanced pairing is referred 
to as "S'C/(3)-pairing" or "super-pairing". [15] We shall 
calculate the gap by writing down the "super" -pair wave 
function. As in single-level pairing, however, we obtain 
the same result by solving the gap equation. 

Assuming SU (3) symmetry the pairing can be calcu- 
lated variationally with a pair wave function that is a gen- 
eralization of Eq. (14) to a sum over I = np,np — 2, 1 
or 



$ (n F ) = 



J2i a i 4>o{nF,l) 

E,«?] 1/2 



(19) 



The weights a; can be found by a variational method 
with the overlap integrals J V? n {R? n , v r 2 dr calculated nu- 
merically. [24] For large np we find a; ~ y/2l + 1 very 
accurately. The sums then reduce to the h.o. level den- 
sity 



£l«r)| 2 = J_£(2Z + l)ft 2 „ 

lm I 



dn 



Here, p a = p/2 is the density of one spin state (see Eq. 
(5) for n = rip) when shells up to energy (n + 3/2)Huj 
are occupied and dp(r)/dn is then the density of the n'th 
shell. 

The sum over 1-states in the gap equation simplifies 
enormously by using the Thomas- Fermi identity valid for 
n > 1 



47T 



dn 



Rrp p 



(20) 



The pairing energy is thus 



g J drr 2 [dp a (r)/dn\ 

r 

2 J drr 2 dpcr(r)/dn[ 



G, 



(21) 



where the supergap G was given in Eq. (9). The con- 
dition for single-shell pairing in atomic h.o. traps is: 
D<2A = 2G, is thus n f <64/5 according to Eqs. (9) 
and (10) or, equivalently, to A^SlO 3 trapped particles. 

For more than two particles in the shell we may again 
invoke the seniority scheme. It applies approximately 
to many supcrpairs when the single-level gap is replaced 
by the supergap G, and the level degeneracy for a full 
h.o. shell £1 = (np + l)(np + 2). Numerical calculations 
for half filled shells confirms that the pairing energy per 
particle is given by the supergap for single-shell pairing 
[16]. 

The single-shell pair energy spectrum does, however, 
differ from that for single-level pairing. In the latter 
case, one can generally write the pairing wave function in 
Eq. (14) as any linear combination of the time-reversed 
two-particle states ipi m (ri)ipi- m (r 2 ). Due to rotational 
SO(3) symmetry the overlap integrals in a single-level 
are independent of m and the spectrum of pairing ener- 
gies is particularly simple: one state with 2E np i as given 
by Eqs. (16-18) and rest of the (21 + 1) states has zero 
energy per pair. In the shell, however, the overlap be- 
tween the radial wave functions depend on I and I' the 
corresponding matrix for I, V — np, np — 2, 1 or has 
a non-trivial eigen-value spectrum: Es — 2G/(2S + 1), 
S = Q,l,...,Int((n F + l)/2), whcnn F > 1 [24]. The cor- 
responding eigen-vectors ai are very well approximated 
by Chebyshev's polynomia as function of x = (l/np) 2 in 
the large np limit. Consequently, the excitation energy 
of a pair is only 2/3 of the energy 2G that it takes to 
break the pair completely. The richer spectrum is im- 
portant for expressing an effective Hamiltonian for BCS 
pairing. 



C. Multi-shell pairing 

For increasing interaction strength pairing also take 
place between different shells, i.e. n ^ np contributes 
in the sum over n in the gap equation (11). The pair- 
ing strength A„ depends only weakly on n for the shells 
around n — np which mainly contribute to the pairing 
and we therefore approximate A„ ~ A„ F . For half filled 
shell p = e nF i and A„ F j = En f i = A. The gap equation 
then becomes [16] 



A = G + G Y - (22) 
„X ^[(n-np)h^ + A 2 

The sum can be converted into an integral when correct- 



ing the lower limit n = 1 by 7 



where C — 0.577.. is 



Eulcr's constant. It then gives to leading orders 



4 



A ^ G + 2G ^ ln 2 ^ 

^ \ 1 + y/l + {"fA/flLo) 2 



which for A <§C ?kj gives 



A = 



G 



l-21n( 7 n F )^ 



(23) 



(24) 



This is the multi-shell pairing gap valid for interactions 
such that Gln(n F )/fiui £1/2. The region of multi-shell 
pairing 0.l£G\n(n F )/huj£,l/2 extends the single-shell 
pairing gap G up to stronger interactions where A <~ Huj 
(see Fig. 1). Cooper pairs are still essentially only formed 
between states within the same shell. However, there is 
pairing in many shells besides that at the Fermi level 
resulting in a gap that is larger than A = G. 



number is replaced by its spatial average in a finite sys- 
tem: 



(k F (r)) 



Jk 2 F (r)d 3 r _ _32_ 
J k F (r)d 3 r I5ir 



MO), 



(27) 



where k F (0) = \J2n F / 1 a osc . That the spatial average 
should be performed over k F (r) 2 follows from Eq. (9) 
where the factor dp(r)/dn oc k F (r) of Eq. (20) en- 
ters twice. By definition of the supergap we find that: 
hw/G = n/2\a\ (k F (r)}, and the two exponents match. 
The prefactors differ due to the approximate cutoff 
n£,2n F in the gap equation, the approximate treatment 
of the of the overlap integrals and their n-dependence in 
the gap equation leading to Eq. (26). 



D. Bulk superfluidity 

For a large system the local supcrfluid field is related 
to that in a uniform system. In the TF approximation 
the local pairing field then depends on radius 

*w "4s*- <25) 

where k F (r) = (37r 2 p(r)) 1/3 = yj2n F - r 2 / a 2 osc / a osc is 
the TF wave number. The prefactor is k = 8/e 2 without 
and k — (2/e) 7 / 3 with corrections from induced interac- 
tions [25]. 

When the attractions are sufficiently strong and/or the 
number of particles in the trap large enough such that 
G\n(-yn F )^hu>, the pairing field exceeds the h.o. fre- 
quency in the bulk of the trap except in a narrow surface 
region. For that reason the low lying single particle exci- 
tations [26,16] and the collective modes are surface modes 
with typical excitation energies of order ~ Two [27]. For 
example, the lowest s-wave excitation (the monopole) is 
uj = 2co in the collisionless limit. This is contrast to the 
weakly interacting system, i.e. below the bulk superflu- 
idity region of Fig.l, where ui a = A for open shells. We 
refer to Ref. [16,27] for a discussion of such modes. 

Here, we concentrate on the bulk supcrfluid field inside 
the trap and will relate it to that in a uniform system. 
Contrary to the surface states, the excitations at higher 
energies are less affected by the pairing field and the 
Bogoliubov wave functions are approximately given by 
the unperturbed h.o. wave functions. Consequently, the 
intra-shell paring condition discussed in Ref. [16], which 
leads to the gap equations (11) and (22), are approx- 
imately valid for these states. However, the gap is to 
be understood as an average superfluid gap (A) in bulk. 
The gap equation (23) is readily solved when (A}^hw 

(A) ~ 2n F hLU cxp(-ftu/2G) . (26) 

We now observe that the exponent in Eq. (26) coincides 
with that in the TF-BCS Eq. (25) when the Fermi wave 



E. Dense liquid 

The dense or strongly interacting limit, fci?|a|^l, can 
be encountered near a Feshbach resonance where a — > 
— oo. In this regime the dilute gas approximation implicit 
in the interactions in the Hamiltonian (1) is no longer 
valid and a new scaling region appears [6-8,12,13]. Both 
the energy per particle and the pairing gap approach a 
finite fraction of the Fermi energy. Recent Greens func- 
tion Monte Carlo calculations [13] find that the inter- 
action energy is -0.56 times the kinetic energy and that 
the odd-even staggering energy or pairing gap in bulk is 
A ~ 0.54£ F . 

F. Multi-level pairing 

We finally address the intermediate regime, G^D^A, 
(see Fig. 1), which lies between the single-level, single- 
shell and multi-shell pairing regions studied above. Here, 
pairing takes place between multiple (I, I') levels lying 
close to the Fermi surface /i ~ £ nF .i, and will therefore be 
referred as multi-level pairing. This region overlaps into 
the regions studied in the above subsections and therefore 
applies to most systems of trapped fermions except traps 
with very few (single-shell), very many (single- level) or 
strongly interacting fermions (multi-shell). As we shall 
see below nuclei can also be considered to belong to the 
multi-level regime. 

For large level splitting D » A the detailed structure 
of the overlap integrals in Eq. (11) becomes important. 
We shall calculate these analytically using WKB as well 
as numerically from the exact h.o. radial wave functions. 
Inserting the latter in the gap equation we calculate the 
pairing gaps as shown in Figs. 2+3 for half filled levels, 
i.e., for fj, = e np i, where E npl = A„ F | = A. 
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FIG. 2. The pairing gap A in units of G vs. the number 
of particles (or uf) in a trap with half filled levels calcu- 
lated from the gap equation (11). The scattering length is 
a — — 0.01a osc . The pairing gaps A nj?i are plotted with + for 
each half filled level I (only even I and hf are shown). The 
gap undergoes a transition from single-shell to multi-level and 
finally single-level pairing. 



A qualitative understanding of these results can be ob- 
tained from analytical WKB calculations. At large uf 
the WKB wave functions TZ np i of Eq. (6) are excellent 
approximations except for the multi-level overlap, I = V 
where a logarithmic singularity appear, which was re- 
sponsible for the ln(i) in Eq. (18). Otherwise, when 
< \l 2 — l' 2 \ <C rip the overlap integral of the WKB wave 
functions of Eq. (6) is dominated by the overlap near the 
turning points. Defining: £ = 2np(r / a osc ) 2 — (r/a osc ) 4 , 
the overlap is 



Jo 



drr 2 nl Fl (r)nl Fl ,(r)*( 



dr 



oc-ln(\x-x'\), (28) 



to leading logarithmic order, where x = (l/np) 2 and 
x' = (/' /np) 2 . This logarithmic dependence of the ma- 
trix element has the interesting consequence that it leads 
to a double log in the gap equation. When the energy fac- 
tor l/E nV ~ (D 2 (x - x 1 ) 2 + A 2 )" 1 / 2 is summed over V 
or, equivalently when np 3> 1, integrated over x', it leads 
to the (usual) factor ln(D/A). The other logarithm ap- 
pears from the matrix element of Eq. (28) which, when 
summed over I' in the gap equation with the energy fac- 
tor, attains the lower cutoff at \x — x'\ ~ A/D. It is 
assumed that the ^-dependence of the gap is sufficiently 
weak within the levels over which pairing take place, i.e. 
A <C D, which is supported by numerical results (see 
Figs. 2+3). The gap equation then reduces to: 

1 ~ a(x)G ^ dx' - - ln \ x - ^l 

Jo \jD 2 \x - x'\ 2 + A 2 



a(x)^ln 2 (p(x)D/A). 



(29) 



to leading logarithmic orders. The factors (3{x) and a(x) 
both depend on x = (l/np) 2 , but not on I or np sepa- 
rately. Only G and D depend on rip explicitly. 
The resulting multi-level gap is 



A = (3(x) D exp \-y/D/Ga(x) 



(30) 



At mid-shell, x = 1/2, the exponent can be calculated: 
a(l/2) ~ 15/32 V2, and the prcfactor: /3(l/2) = 2e, to 
leading logarithmic order for large np. Near open (I = 0) 
and closed (/ = np) shells there are generally only half 
as many states to pair with which reduces a(x — 0) and 
a(x = 1) by a factor ~ 1/2 on average. However, as 
the matrix elements are larger for / ~ np but smaller 
for I = (sec Eq.(28)), the gaps become asymmetric 
with a maximum above midshcll. Such a shell structure 
is also found in the multi-level pairing regime when the 
gap equation is solved numerically with the exact h.o. 
wave functions as is shown in Figs. (2+3). 

Two illustrative examples are shown in Figs. (2+3). In 
the first the coupling is sufficiently weak, a = — 0.01a osc , 
that the pairing undergoes transitions from from single- 
shell at small np to multi-level for lO^np^AO and finally 
single- level pairing for large np. In Fig. 3 the coupling 
is stronger, a = — 0.04a osc , and the pairing undergoes 
transitions from from single-shell to multi-level and ap- 
proaches multi-shell pairing for large np. These transi- 
tions between pairing regimes are illustrated in Fig. 1. 




(3N) 1/3 = n F 

FIG. 3. Same as Fig. 2. but with a = -0.04a osc . Full line 
is the multi-shell gap G/[l — 2 hi(jnp)G /hu>]. The gap under- 
goes a transition from single-shell to a mixture of multi-level 
and multi-shell pairing. 



The multi-level gap of Eq. (30) is quite robust and 
applies to many systems as long as G^D^A. Its validity 
for nuclei will be discussed below. The formula for the 
multi-level gap can be generalized by relating the shell 
splitting to the density of states at the Fermi surface as 
dn/de = Q/D, where fi = np is the number of states 
in the shell. Likewise, the supergap is related to the 
coupling constant as G cx gfl. If the level spectrum of 
Eq. (8) is changed, the multi-level gap remains valid, if 



G 



the shell splitting is correspondingly scaled with the level 
density. In other words the double log does not depend 
on the details of the level spectrum but is generic for 
systems like h.o. traps, because the one log is associated 
with the overlap matrix elements between nearby states 
and the other is associated with the quasi-particle energy 
in the gap equation. 

The appearance of a double log is not unique for multi- 
level pairing but also occurs in the case of color supercon- 
ductivity in quark matter [28]. The physical mechanism 
behind is, however, different. Within pcrturbativc QCD 
the singular quark-quark interaction gQCD/q 2 , where q 
is the momentum transfer carried by a gluon, is dynam- 
ically screened by Landau damping. When the interac- 
tion is integrated over momentum and energy transfer, 
the dynamical screening leads to a logarithm of the gap. 
Thus it is the interactions that are responsible for the 
second logarithm and not the wave function overlap of 
nearby states as in the h.o. trap. 

IV. PAIRING IN NUCLEI 

The nuclear mean field is often approximated by a sim- 
ple h.o. form and the residual effective pairing interac- 
tion by a delta force in order to obtain some qualitative 
insight into single particle levels, pairing, collective mo- 
tion, etc., (see, e.g., [19,23]). We can therefore compare 
pairing in nuclei to that in traps as investigated above, 
once the h.o. potential is adjusted to describe nuclei. We 
emphasize that we do not intend to calculate the quan- 
titative pairing gaps for each individual nucleus which 
would require detailed knowledge of the individual level 
spectra, deformation, many-body effects, etc. Instead we 
aim at qualitative results for the pairing gap dependence 
on mass number, shell effects, and to extrapolate to very 
large nuclei, nuclear and neutron star matter. 

Large nuclei have approximately constant central den- 
sity po ~ 0.14 fm~ 3 and Fermi energy Ep in bulk. 
Therefore the h.o. frequency, which is fitted to the nu- 
clear mean field, decreases with the number of nucleons 
A = N + Z, where N now is the number of neutrons and 
Z the number of protons in the nucleus, as 

hu ~ Ep/np ~ 41MeV • A~ 1/3 . (31) 

In the valley of /^-stability the number of protons is Z ~ 
A/(2 + 0.0155A 2 / 3 ). 

Secondly, the nuclear mean field deviates from a h.o. 
potential by being almost constant inside the nucleus and 
vanish outside. The resulting net anharmonic nuclear 
field is stronger and opposite in sign to the corresponding 
(anharmonic) mean field in atomic traps. Therefore, the 
level splitting is larger and the ordering of the 1-levels 
is reversed. In addition, a strong spin-orbit force splits 
the single particle states of total angular momentum j = 
I ± 1/2, such that the j = rip + 1/2 is lowered down to 
the shell (rip — 1) below. 



Proton and neutron pairing gaps are typical of order 
~ 1 MeV in nuclei, i.e., smaller than both Tiui and D 
but of order the average splitting between two adjacent 
j-levels. Consequently, nuclei can be considered as h.o. 
traps with a level splitting D such that they fall into the 
multi-level pairing regime. However, the number of par- 
ticles in nuclei is relatively small that nuclei are close to 
the single shell pairing regime. Furthermore, the interac- 
tions are so strong that multi-shell pairing also becomes 
important. From the results of Sec. Ill we can predict 
several features of neutron and proton pairing gaps: 

• Mass scaling: Since Huj scale as ~ A -1 / 3 and 
a osc oc rip 2 (see Eq. (9)) the single-shell pairing 
gap also scales as G ~ A -1 / 3 . The level split- 
ting tend to reduce the pairing towards the multi- 
level gap but is compensated by multi-shell pairing. 
Therefore, the pairing gaps in light and medium 
mass nuclei scale approximately as: A ~ G cx 
A-V3. 

• Shell structure: The pairing gaps should ex- 
hibit a strong shell structure similar to those for 
multi-level pairing (Figs. 2+3), however with re- 
versed I since D is negative for nuclei. Due to 
the strong spin-orbit force the j = I ± 1/2 states 
are split and the j = rip + 1/2 is lowered down 
to the shell below. The magic numbers become 
N, Z = 8, 14, 28, 50, 82, 126, 184, .., etc. rather than 
the h.o. filled shell particle numbers N,Z = 
2, 8, 20, 40, 70, 112, 168, 240, etc. 

• Bulk limit for large nuclei/nuclear matter: 

For very large nuclei multi-shell pairing becomes 
important and pairing approaches that in bulk mat- 
ter. By fitting the effective nucleon coupling con- 
stant to pairing gaps in finite nuclei we will be- 
low estimate the pairing in nuclear matter from Eq. 
(25). 

These predictions agree qualitatively with experimental 
data (see Figs. 4+5). 

For a more quantitative calculation the level spectrum 
must be specified. Instead of fitting the level spectrum 
for each nucleus with a correspondingly large number of 
adjustable parameters, we make the following simplifying 
approximation analogous to the single particle spectrum 
ofEq. (8) 

The level splitting is approximated by D — 0.13(np + 
3/2)tku in analogy with Eq. (10) and based on nuclear 
structure calculations, which generally find that the level 
splitting increase with shell number and is of order D ~ 
Tiw for heavy nuclei (np ~ 6). The spin-orbit splitting 
and resulting change in magic numbers are incorporated 
to first approximation by simply moving the lowest level 
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in each shell down to the shell below, i.e., D = Tiuj for 
j = n F + 1/2. 

The pairing gaps and quasi-particle energies can now 
be calculated by solving the gap equation inserting the 
h.o. matrix elements and with Jkv and level splitting as 
given by Eqs. (31) and (32). The effective strength a is 
then the only adjustable parameter. 

The data on neutron and proton pairing is obtained 
from the odd-even staggering of nuclear binding energies 
B(N, Z). It has been shown that mean field contributions 
can be removed [30] by using the three-point filter 

<-l) N 

A (3) (A0 = v ; [B(N -1,Z) + B(N + 1, Z) 

-2B(N,Z)}, (33) 

when N is an odd number of neutrons. The analogous 
relation is valid for protons. 

The total binding energies are elaborate sums over 
quasi-particle energies weighted with occupation num- 
bers. Mean field energies and single particle energies 
should also be included self consistently. Furthermore, 
the Bogoliubov transformation does not preserve exact 
particle number and does not treat a single un-paired par- 
ticle that appears in odd-number systems. Also nuclear 
spectra are complicated by deformations and the finite 
range, spin and spin-orbit dependences of the nucleon- 
nucleon interaction [19,20]. For these reasons we can- 
not calculate A' 3 ) directly. Instead we will for simplicity 
compare to the quasi-particle energies and pairing gaps, 
which are calculated directly from the gap equation as 
function of fi(N,Z). It has been argued (see, e.g. [30] 
and refs. herein) that mean-field effects cancel in the 
quasi-particle energy for odd particle number and that it 
therefore may be compared to the corresponding three- 
point energies. However, as particle number fluctuations, 
deformations, possible mean field energy corrections and 
other effects are not included, the pairing gaps are also 
shown in Figs. (4+5) for comparison. They are equal 
to the quasi-particle energy for half-filled levels only and 
generally smaller especially at the magic numbers where 
the gaps vanish in several cases. 

We compare in Fig. 4 to the experimental A^ 3 ^(7V) 
averaged over isotopes with the calculated gaps A„ F ;, 
and in Fig. 5 the analogous for protons A^ 3 \Z) averaged 
over isotones. In the calculations the effective coupling 
is the only parameter fitted to experimental data. For 
both neutrons and protons we extract a ~ —0.41 fm. 




FIG. 4. Neutron pairing energies vs. the number of 
neutrons. The experimental odd-even staggering energies 
A^ 3 \N) are averaged over isotopes [31,30]. The calculated 
gaps A and quasi-particle energies E are obtained from 
the gap equation (see text) with effective coupling strength 
a — —0.41 fm. The supergap G is shown with dashed line. 

We note that although the neutron pairing gaps in Fig. 
4 generally are larger than the proton ones in Fig. 5, this 
is not reflected in the effective coupling constants. The 
reason is the asymmetry of heavy nuclei. For example, 
for N — 82 the mass number is A ~ 140 whereas for 
Z = 82 it is A ~ 208. The mass numbers enter both 
D and G and lead to a reduction of the proton pairing 
gap relative to the neutron one by just the right amount 
so that the experimental data on neutron and proton 
odd-even staggering can be fitted with the same pairing 
strength a — —0.41 fm. 

Considering the simplicity of the model it describes a 
large number of experimental gaps fairly well on average. 
In a number of cases, however, the calculated pairing 
gaps differ significantly from the measured neutron gaps. 
Some of these deviations can be attributed to the crude 
single particle level spectra assumed. If the single particle 
level energies are adjusted according to more detailed 
mean field calculations (see, e.g., [19]) the agreement with 
experimental pairing gaps improves in several cases. 
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FIG. 5. As Fig. 4 but for protons 



The pairing gaps are sensitive to the shell splitting and 
the coupling. The uncertainty in the coupling is smaller 
because the changes in D affects both the exponent and 
the prefactor in the gap of Eq. (30) in a compensating 
way. Another uncertainty arise from the upper cutoff. 
Whereas h.o. traps can pair between shells up to ~ 2np, 
nuclei have a continuum of states at about the binding 
energy per nucleon Eb — 8 MeV above the Fermi level. 
Therefore the sum over shells in the gap equation has 
been limited to n^(Ep + Eb)/?iu>, which for medium 
mass nuclei corresponds to n^riF + 1, with a smooth 
cutoff as in Ref. [32]. The pairing gaps are, however, 
only logarithmically sensitive to this cutoff. 

Nuclei can also be deformed around mid-shell which in- 
crease pairing [30]. These and other nuclear many-body 
effects must be included in a more quantitative study of 
pairing in individual nuclei. For this purpose it will be 
useful also to study pairing in elongated atomic traps 
with deformations like in nuclei. Near closed shells pair- 
ing reduce or inhibit deformations. The effect of pairing 
is also observed in rotational spectra where the moment 
of inertia is reduced from the rigid to superfluid or irro- 
tational value [19]. In traps the external h.o. potential 
generally dominates over the mean field interaction en- 
ergy and thus does not deform spontaneously [33] as in 
nuclei - unless of course if the h.o. trap itself is deformed. 

Nuclei cannot directly be placed in any of the various 
pairing phases of Fig. 1 because the level splitting is 
larger than for trapped atoms. The effective scattering 
length a = —0.41 fm and a osc ~ lfm/y/nF would place 
nuclei with masses up to A = 250 corresponding to rip ,6 6 
in the upper left corner of Fig. 1 in the multi-shell re- 
gion. However, because D is much larger in nuclei than 
for trapped atoms, the multi-level pairing region extends 
down to lower uf and up to larger strengths \a\. Fur- 
thermore, the continuum of states in nuclei reduce the 
effect of multi-shell pairing as discussed above. There- 
fore, nuclei rather belong to the transition region between 
single-shell, multi-shell and multi-level pairing. 

The best fit to odd-even staggering energies of nuclei 
determines a = —0.41 fm accurately. Systematic errors 
may, however, be expected from the approximations im- 
plicit in the level splitting, the cut-off, and in approximat- 
ing A' 3 ) by A. However, because the multi-level and the 
multi-shell pairing partly compensate, a good approxi- 
mation to the average gap in nuclei is the single-shell 
supergap 



A ~ G 



\a\ 5.5MeV 
0.41fm A 1 / 3 



(34) 



This supergap is also shown in Figs. (4+5). It docs not 
depend on the level-splitting or cut-off and is therefore 
a robust prediction for the average magnitude and mass 
scaling of pairing gaps in nuclei. 



Empirically the pairing term in Bethe-Weiszacker 
liquid-drop formula, A ~ 12 MeV-A -1 / 2 , fits the odd- 
even staggering energies of nuclei with A^,250 averaged 
over shell effects. The scaling with mass number can now 
be understood in terms of the supergap with shell correc- 
tions. If D was a constant times Tiuj then the multi-level 
gap of Eq. (30) would also scale as A -1 / 3 . However, be- 
cause D increase with np oc A -1 / 3 , the multi-level gap 
decrease faster with A for small and medium mass nu- 
clei but slower for heavy nuclei due to multi-shell pairing. 
Both are in accordance with the empirical mass depen- 
dence. 

The pairing in nuclear matter can also be estimated 
once the effective interaction has been determined. In- 
serting in Eq. (25) a = —0.41 fm and Uf = 1-33 fm -1 
at nuclear saturation density, p = 0.15 fm -1 , we obtain 
the proton and neutron pairing gaps 



A ~ l.lMeV, 



(35) 



in the bulk of very large nuclei and in symmetric nu- 
clear matter at nuclear saturation density. This num- 
ber is compatible with earlier calculations [34] of the 1 So 
pairing gap in nuclear and neutron star matter around 
normal nuclear matter densities. A smaller value for bulk 
pairing might have been expected from Figs. (4+5) by 
extrapolating the odd-even staggering energies of heavy 
nuclei to higher mass numbers. However, as A oc n F be- 
comes large the multi-shell pairing contribute with the 
increasing term log(riF) in the gap equation. Therefore 
the bulk value is larger than the pairing gap in heavy 
nuclei which also have a smaller cutoff due to continuum 
states as discussed above. 

Neutron star matter has a wide range of densities and 
is very asymmetric, Z/A ~ 0.1. One can attempt to esti- 
mate of the pairing gaps as function of density from the 
gap in bulk, Eq. (11), with a ~ —0.41 fm and the neu- 
tron or proton Fermi wave numbers kp' Z — (3ir 2 pn.z) 1 ^ 3 
as function of densities. However, the effective interac- 
tion a is density dependent. At higher densities we ex- 
pect the effective interaction to become repulsive as is 
the case for the nuclear mean field at a few times nuclear 
saturation density. At lower densities the effective scat- 
tering length should approach that in vacuum which for 
neutron-neutron scattering is a( 1 S'o) — —18 fm. This di- 
lute limit fejr|a| <C 1 does, however, require extremely low 
densities as compared to normal nuclear matter density. 



V. SUMMARY 

Pairing gaps have been calculated for ultracold atomic 
Fermi gases in harmonic oscillator traps and compared to 
nuclei. The pairing mechanism was found to be similar 
for these systems in the sense that the spacing between 
single particle states (and shells) reduce the pairing over 
several of these levels near the Fermi surface referred to 
as multi-level pairing. At low particle densities the shell 
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structures in traps are pronounced as they are in nuclei 
and the level degeneracies are important for the size of 
the gaps which can differ substantially from those known 
from homogeneous systems [25] and systems with contin- 
uous level densities. 

Neutron and proton pairing gaps in nuclei were cal- 
culated and with an effective coupling strength a — 
—0.41 fm a qualitative description of their shell struc- 
ture could be given, and the average pairing gaps were 
found to scale with mass number approximately as A ~ 
5.5 McV/A 1 / 3 as predicted from the supergap. Eventu- 
ally for large mass number the gap approaches the 1 So 
superfluid gap in uniform nuclear matter which was cal- 
culated from the as A ~ 1.1 MeV for both neutrons and 
protons. 

Mixing fermionic with bosonic atoms improve cooling 
[36,9] to lower temperatures so that weak pairing can 
also be studied, and the additional induced interactions 
between fermions and bosons generally enhance pairing 
[25]. Furthermore, the shell splitting can be changed in 
a controlled way by the number of bosons in the trap 
and the sign and strength of their interaction with the 
fermions. 

The similarity of multi- level and bulk pairing in atomic 
traps and nuclei may provide new insight into pairing 
and superfluidity in nuclei and neutron star matter from 
tabletop experiments at low temperatures. 



In the multi-level regime T c can be determined from 
Eq. (Al) with overlap integrals as given in Eq. (28). To 
leading logarithmic order we find for large uf 

T c = — A(T = 0) . (A4) 

TT 

That T C /A(T = 0) = j/tt as in the uniform Fermi gas is 
mainly because the pairing take place over several 1-levels 
and the level density therefore is effectively continuous. 
The overlap integrals do not change this ratio to leading 
logarithmic order. 

Although T C /A(T = 0) = 7/tt ~ 0.567 is close numer- 
ically to the value 1/2 found in the single-shell, single- 
level and multi-shell regimes, the difference reveals the 
qualitative differences in the underlying level spectrum, 
namely, continuous vs. discrete respectively. 

Near a Feshbach resonance the strongly interacting 
Fermi gas becomes unstable towards molecule forma- 
tion and T c for BCS superfluidity is expected cross- 
over towards the slightly smaller critical temperature for 
forming a Bose-Einstein condensate of molecules [37]. 
Both the BCS and BEC critical temperatures are, how- 
ever, above the lowest temperatures achieved recently for 
trapped Fermi atoms [5,6,9] if we assume T c ~ 0.5A as 
in Eqs. (A2-A4) and take A = 0.54E F according to Ref. 
[13]. 



APPENDIX A: CRITICAL TEMPERATURES 

The pairing gaps generally decrease with increasing 
temperature from its zero temperature value A(T = 0), 
that was calculated above, to the critical temperature T c , 
where it vanishes. The temperature dependence of the 
gap and T c itself are determined from 

A„w(T) = £[1 - 2f{E n ,, v /T))^^lj- 

n I 

/' oc 

x / drrXjto^'M, (Al) 
Jo 

with E n , v = y/(e nn > -/i) 2 + A 2 (T). 

As shown in [29] T c is exactly half of the zero temper- 
ature gap 

T c = G/2(l - 2 ln( 7 n F )G/M = ^A(T = 0), (A2) 

in the single-shell and multi-shell pairing regimes, and 
this also applies to the single-level pairing regime. 

In uniform Fermi gases the critical temperature is [25] 

T c = -kE f exp(2/7rafc F ) = -A(T = 0) . (A3) 

TT TT 

The ratio T C /A(T = 0) = 7/tt ~ 0.567 is the same irre- 
spective of whether induced interactions are included or 
not. 
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